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We present first results for the masses of positive and negative parity excited baryons calculated 
in lattice QCD using an 0(a 2 )-improved gluon action and a fat-link irrelevant clover (FLIC) fermion 
action in which only the irrelevant operators are constructed with APE-smeared links. The results 
are in agreement with earlier calculations of N* resonances using improved actions and exhibit a clear 
mass splitting between the nucleon and its chiral partner. An correlation matrix analysis reveals 



two low-lying J 



states with a small mass splitting. The study of different A interpolating 



fields suggests a similar splitting between the lowest two Ai octet states. However, the empirical 
mass suppression of the A* (1405) is not evident in these quenched QCD simulations, suggesting 
a potentially important role for the meson cloud of the A* (1405) and/or a need for more exotic 
interpolating fields. 
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I. INTRODUCTION 

Understanding the dynamics responsible for baryon ex- 
citations provides valuable insight into the forces which 
confine quarks inside baryons and into the nature of QCD 
in the nonperturbative regime. This is a driving force 
behind the experimental effort of the CLAS Collabora- 
tion at Jefferson Lab, which is currently accumulating 
data of unprecedented quality and quantity on various 
N — > N* transitions. With the increased precision of 
the data comes a growing need to understand the ob- 
served N* spectrum within QCD. Although phenomeno- 
logical low-energy models of QCD have been successful 
in describing many features of the N* spectrum (for a 
recent review see Ref. [lj), they leave many questions 
unanswered, and calculations of N* properties from first 
principles are indispensable. 

One of the long-standing puzzles in spectroscopy has 
been the low mass of the first positive parity excitation 
of the nucleon (the J p = i + A r *(1440) Roper resonance) 
compared with the lowest lying odd parity excitation. In 
a valence quark model, in a harmonic oscillator basis, 



the \ 



i + 



state 



state naturally occurs below the N — 2, 
0. Without fine tuning of parameters, valence quark 
models tend to leave the mass of the Roper resonance 
too high. Similar difficulties in the level orderings ap- 
pear for the | + A*(1600) and ± + S*(1690), which has 
led to speculations that the Roper resonances may be 
more appropriately viewed as "breathing modes" of the 
states |3j , or described in terms of meson-baryon dynam- 
ics alone 0, or as hybrid baryon states with explicitly 
excited glue field configurations Q . 

Another challenge for spectroscopy is presented by the 
A 1 / 2- (1405), whose anomalously small mass has been in- 



terpreted as an indication of strong coupled channel ef- 
fects involving T,tt, KN, ... 0, and a weak overlap with 
a three-valence constituent-quark state. In fact, the role 
played by Goldstonc bosons in baryon spectroscopy has 
received considerable attention recently 0, Q • 

It has been argued 9] that a spin-flavour interaction 
associated with the exchange of a pseudoscalar nonet of 
Goldstone bosons between quarks can better explain the 
level orderings and hyperfine mass splittings than the 
traditional (colour-magnetic) one gluon exchange mech- 
anism. On the other hand, some elements of this ap- 
proach, such as the generalisation to the meson sector 
or consistency with the chiral properties of QCD, re- 
main controversial |llllOllll|. Furthermore, neither spin- 
flavour nor colour-magnetic interactions are able to ac- 
count for the mass splitting between the A 1 / 2 (1405) 
and the A 3 / 2 (1520) (a splitting between these can arise 
in constituent quark models with a spin-orbit interac- 
tion, however, this is known to lead to s pur ious mass 
splittings elsewhere pi H^|L Recent work |l3j on nega- 
tive parity baryon spectroscopy in the large- N c limit has 
identified important operators associated with spin-spin, 
spin-flavour and other interactions which go beyond the 
simple constituent quark model, as anticipated by early 
QCD sum-rule analyses |l4j . 

The large number of states predicted by the con- 
stituent quark model and its generalisations which have 
not been observed (the so-called "missing" resonances) 
presents another problem for spectroscopy. If these states 
do not exist, this may suggest that perhaps a quark- 
diquark picture (with fewer degrees of freedom) could 
afford a more efficient description, although lattice simu- 
lation results provide no evidence for diquark clustering 
[l5| . On the other hand, the missing states could sim- 
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ply have weak couplings to the ttN system pj. Such a 
situation would present lattice QCD with a unique op- 
portunity to complement experimental searches for N*'s, 
by identifying excited states not easily accessible to ex- 
periment (as in the case of glueballs or hybrids). 

In attempting to answer these questions, one fact that 
will be clear is that it is not sufficient to look only at the 
standard low mass hadrons (it, p, N and A) on the lattice 
— one must consider the entire N* (and in fact the entire 
excited baryon) spectrum. In this paper we present the 
first results of octet baryon mass simulations using an 
0(a 2 ) improved gluon action and an improved Fat Link 
Irrelevant Clover (FLIC) ^| quark action in which only 
the irrelevant operators are constructed using fat links 
|17|. Configurations are generated on the Orion super- 
computer at the University of Adelaide. After reviewing 
in Section II the main elements of lattice calculations of 
excited hadron masses and a brief overview of earlier cal- 
culations, we describe in Section III various features of 
interpolating fields used in this analysis. Section IV re- 
views the details of the lattice simulations, and Section V 
gives an overview of the methodology for isolating baryon 
resonance properties. In Section VI we present the results 
from our simulations and in Section VII make conclud- 
ing remarks and discuss possible future extensions of this 
work. 



lowest octet and decuplet positive and negative parity 
baryons using overlap fermions with pion masses down 
to ~ 180 MeV. While these authors claim to have ob- 
served the Roper in quenched QCD, it remains to be 
demonstrated that this conclusion is independent of the 
Bayesian-prior assumed in their analysis lisl l2fil |. 

Following standard notation, we define a two-point cor- 
relation function for a spin-^ baryon B as 

G B (t,p) = e~** <n \xb(x) X b(0)\ fi) (1) 

where \b is a baryon interpolating field and where we 
have suppressed Dirac indices. All formalism for corre- 
lation functions and interpolating fields presented in this 
paper is carried out using the Dirac representation of 
the 7-matrices. The choice of interpolating field \b is 
discussed in Section III below. The overlap of the inter- 
polating field xb with positive or negative parity states 
\B ± ) is parameterised by a coupling strength X B ± which 
is complex in general and which is defined by 



<fi | XB (0) \B+,p, s) = A B+ y u B+ (p, s) , (2a) 



(SI\ X b(0)\B-, P ,s) = \ b -, 



Er 



lhU B -(p,s) ,(2b) 



II. EXCITED BARYONS ON THE LATTICE 

The history of excited baryons on the lattice is quite 
brief, although recently there has been growing interest 
in finding new techniques to isolate excited baryons, mo- 
tivated partly by the experimental N* program at Jef- 
ferson Lab. The first detailed analysis of the positive 
parity excitation of the nucleon was performed by Lein- 
weber |l8[ using Wilson fermions and an operator prod- 
uct expansion spectral ansatz. DeGrand and Hecht 0] 
used a wave function ansatz to access P-wave baryons, 
with Wilson fermions and relatively heavy quarks. Sub- 
sequently, Lee and Leinweber |20j introduced a parity 
projection technique to study the negative parity \ 
states using an C(a 2 ) tree-level tadpole- improved D X 34 
quark action, and an 0(a 2 ) tree-level tadpole-improved 
gauge action. Following this, Lee [2lJ reported results us- 
ing a D234 quark action with an improved gauge action 

1 + 1 — 

on an anisotropic lattice to study the 5 and ^ exci- 
tations of the nucleon. The RIKEN-BNL group [13 has 

also performed an analysis of the N* ( i ) and N' ( i + ) ex- 
cited states using domain wall fermions. More recently, a 
nonperturbatively improved clover quark action has been 
used by Richards et al. [2^ to study the N*(h ) and 
A* (| ) states, while Nakajima et al. have studied the 

7V*(i ) and A*(i ) states using an anisotropic lattice 
with an O(a) improved quark action [24J. Constrained- 
fitting methods based on Bayesian priors have also re- 
cently been used by Lee et al. [25j to study the two 



where M B ± is the mass of the state B ± , E B ± = 

l^f is its energy, and u B± (p, .) is a Dirac spinor 

with normalisation u B± (p, s)u g ±(p, s) = S af3 . For large 
Euclidean time, the correlation function can be written 
as a sum of the lowest energy positive and negative parity 
contributions 



G B (t,p) « A|. 



( 7 -p + M s+ ) 



2E 



A2 (7-P-M b -) 
B 2E B - 



(3) 



when a fixed boundary condition in the time direction is 
used to remove backward propagating states. The posi- 
tive and negative parity states are isolated by taking the 
trace of G B with the operator T + and T_ respectively, 
where 



r± = 



1 + M ^ 
1 ± — — 74 



E 



B I 



(4) 



For p = 0, Tj_ = r± so that T± are then parity projec- 
tors. For p = 0, the energy E B ± = M B ± and using the 
operator T± we can isolate the mass of the baryon B ± . 
In this case, positive parity states propagate in the (1, 1) 
and (2, 2) elements of the Dirac matrix of Eq. J3J|, while 
negative parity states propagate in the (3, 3) and (4, 4) 
elements. 

In terms of the correlation function G B , the baryon 
effective mass function is defined by 

M B (t) = log[G B {t, 0)] - log[G fl (f + 1,0)]. (5) 
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Meson masses are determined via analogous standard 
procedures. 

III. INTERPOLATING FIELDS 

In this analysis we consider two types of interpolating 
fields which have been used in the literature. The nota- 
tion adopted is similar to that of Ref. [23 ■ To access the 
positive parity proton we use as interpolating fields 

Xi + (x) = e a bc (Ua(x) C75 d b (x)) u c (x) , (6) 

and 

Xl + ( x ) = e abc(ul(x) C d b (x))j 5 u c (x) , (7) 



where the fields u, d are evaluated at Euclidean space- 
time point x, C is the charge conjugation matrix, a, b 
and c are colour labels, and the superscript T denotes 
the transpose. These interpolating fields transform as 
a spinor under a parity transformation. That is, if the 
quark fields q a {x) (q = u,d,- ■ • ) transform as 

Vq a {x)V^ = +7o<? a (i) , 
where x — (xq, —x), then 

T X p+ (x)V^ =+ 7 oX p+ (i) ■ 
For convenience, we introduce the shorthand notation 



g(S fl ,S h ,S h ) = e abc e a ' b ' c '\S a f f(x,0)t r S b / 2 T (x, 0)S% (x, 0) + S% (x, 0) S% T (x, 0) S% (x, 0) 



bb' T / rwocc'i 



?bb' T 1 



(8) 



where S , ? a _ 3 (x, 0) are the quark propagators in the background link-field configuration U corresponding to flavours 
/i_3. This allows us to express the correlation functions in a compact form. The associated correlation function for 
Xi + can be written as 



Git{t,P)T) = /j^e-^tr [-r g (S u , CS d G-\ S„ 



(9) 



where (■ • • ) is the ensemble average over the link fields, T is the r± projection operator from Eq. |0J, and C — C75. For 
ease of notation, we will drop the angled brackets, (•••), and all the following correlation functions will be understood 
to be ensemble averages. For the \2 + interpolating field, one can similarly write 



G p +(t,P;T) =^e-^tr [-r Q (75^75, CS d C~\ 75^75 

while the interference terms from these two interpolating fields are given by, e.g., 

G p +(t,p-T) - ^e-^tr -T U (s u7s , CS d C~\ S ul5 

X 

G p +(t,p-T) = ^ e - l ^tr[-r \g(-y 5 s u , CS d C~\ l5 S u 



(10) 

(11) 
(12) 



The neutron interpolating field is obtained via the ex- 
change u «-> d, and the strangeness -2, 5 interpolating 
field by replacing the doubly represented u or d quark 
fields in Eqs. © and J7J) by s quark fields. S and S 
interpolators are discussed in detail below. 

As pointed out in Ref. , because of the Dirac struc- 
ture of the "diquark" in the parentheses in Eq. ©, in the 
Dirac representation the field Xi + involves both products 
of upper x upper x upper and lower x lower x upper com- 
ponents of spinors for positive parity baryons, so that in 
the nonrelativistic limit Xi + — C(l)- Here upper and 



lower refer to the large and small spinor components in 
the standard Dirac representation of the 7 matrices. Fur- 
thermore, since the "diquark" couples to total spin 0, one 
expects an attractive force between the two quarks, and 
hence better overlap with a lower energy state than with 
a state in which two quarks do not couple to spin 0. 

The x P+ interpolating field, on the other hand, is 
known to have little overlap with the nucleon ground 
state 0, |28| . Inspection of the structure of the Dirac 
matrices in Eq. Q reveals that it involves only prod- 
ucts of upper x lower x lower components for positive 
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parity baryons, so that X2 + = 0{p 2 /E 2 ) vanishes in the 
nonrclativistic limit. As a result of the mixing of up- 
per and lower components, the "diquark" term contains 
a factor a ■ p, meaning that the quarks no longer couple 
to spin 0, but are in a relative L = 1 state. One ex- 
pects therefore that two-point correlation functions con- 
structed from the interpolating field xf + are dominated 
by larger mass states than those arising from Xi + at early 
Euclidean times. 

While the masses of negative parity baryons are ob- 
tained directly from the (positive parity) interpolating 
fields in Eqs. © and Q by using the parity projec- 
tors T±, it is instructive nevertheless to examine the 
general properties of the negative parity interpolating 
fields. Interpolating fields with strong overlap with the 
negative parity proton can be constructed by multiply- 
ing the previous positive parity interpolating fields by 75, 
Xx~2 = 75 Xi~2- I n contrast to the positive parity case, 
both the interpolating fields Xi~ an d X2~ mrx upper and 
lower components, and consequently both \i~ an d X2~ 
are 0{p/E). 

Physically, two nearby J p = | states are observed 
in the excited nucleon spectrum. In simple quark mod- 
els, the splitting of these two orthogonal states is largely 
attributed to the extent to which scalar diquark config- 
urations compose the wave function. It is reasonable to 
expect Xi to have better overlap with scalar diquark 
dominated states, and thus provide a lower effective mass 



in the moderately large Euclidean time regime explored 
in lattice simulations. If the effective mass associated 
with the X2~ correlator is larger, then this would be evi- 
dence of significant overlap of X2~ with the higher lying 
states. In this event, a correlation matrix analysis 
(see Section V) will be used to isolate these two states. 

Interpolating fields for the other members of the 
flavour SU(3) octet are constructed along similar lines. 
For the positive parity E° hyperon one uses 



Xf(x) = ~^ e abc { (Ua(x) C75 S b (xj) d c (x) 

+ (d T a {x) C75 s b (x))u c (x)} ,(13) 

Xf(x) = -^tabc { (Ua (x) C S b (x)) 75 d c (x) 

+ (dl(x)Cs b (x)) l5 u c (x)} .(14) 

Interpolating fields used for accessing other charge states 
of S are obtained by d — > u or u — > d, producing corre- 
lation functions analogous to those in Eqs. © through 
ifTTJl. Note that xf transforms as a triplet under SU(2) 
isospin. An SU(2) singlet interpolating field can be con- 
structed by replacing " + " — > " — " in Eqs. iTH3|) and 
(THIIl . For the SU(3) octet A interpolating field (denoted 
by "A 8 "), one has 



Xi S {x) = -^=e afcc {2 (ua(x) C75 d b {x)) s c {x) + (ul(x) Cj 5 s b (x)) d c (x) - (d r (x) C75 s b (x)) u c {x)} , (15) 
X2 8 ( x ) = ~^ 6abc { 2 ( u a( x ) Cd b {x)) 75 s c (x) + (ul (x) Cs b (xj) 75 d c (x) - (dl (x) C s b {x)) 75 u c (x) j , (16) 



which leads to the correlation function 



(17) 



G^(t,p;T) = i^V^-tr -r j 2Q (s„ CS u C-\ S d ) +2£ CS d C-\ S u ) 

X 

+ 2Q (^Sd, CSuC" 1 , S^j + 2Q (^S u , CSdC^ 1 , S t 
— G (^Sd, CS S C 1 , sSj — Q {^Su, CS S C 1 , S c 

and similarly for the correlation functions G22 , G\2 an d G21 ■ 

The interpolating field for the SU(3) flavour singlet (denoted by "A 1 ") is given by [27j 

Xi\ x ) = - 2 e afcc| - («„ (x) C75 d b (x)) s c (x) + (v%(x) C75 s b (x)) d c (x) - (d^(x) C75 s b (x)) u c (x)} , (18) 
X2 (x) = -2 e a6c | - (v%(x) C d b (x)) 75 s c (x) + (u^(x) C s b (x)) 75 d c (x) - {d T a (x) C s b (x)) 75 u c (x)} , (19) 



where the last two terms are common to both Xi an d Xi ■ The correlation function resulting from this field involves 
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quite a few terms, 



+ 75 S7 CS C U C *C- 1 Sf 75 + 75^ QQ CS? r C~' Sf 75 + TsST CS S CC ' T C~ 1 S bb '75 
+ 75^ a ' CSfrd- 1 75 ~ 7 5 5r' 7 5 tr 



Cf66' s~i qcc' T 



75 SS a 75 tr 



75^ a 75 tr 



Cjbb' £j qcc' T £j — 1 



(20) 



In order to test the extent to which SU(3) flavour symmetry is valid in the baryon spectrum, one can construct 
another A interpolating field composed of the terms common to A 1 and A 8 , which does not make any assumptions 
about the SU(3) flavour symmetry properties of A. We define 



Xi\x) = -^=e abc { (v%(x) C7 5 s b (x)) d c (x) - (d^(x) C75 s b (x)) u c (x)j , 
x£ C ( x ) = ^ f afc{ ("IWCsbW) 75d c (x) - (dl (x)Csb(x)) 7 5 u c (x)| , 



(21) 
(22) 



to be our "common" interpolating fields which are the isosinglet analog of xf an d X2 m ^qs. (|13[) and l|14|). Such 
interpolating fields may be useful in determining the nature of the A* (1405) resonance, as they allow for mixing 
between singlet and octet states induced by SU(3) flavour symmetry breaking. To appreciate the structure of the 
"common" correlation function, one can introduce the function 



G(S fl , S h ,S h ) = e abc e a '»' c ' Sf(x, 0) tr S b / 2 T (x, 0)S% (x, 0) - Sff (x, 0) S% T (x, 0) S% (x, 0) 



jfcb'T 



bb'Tl 



(23) 



which is recognised as Q in Eq. (JHJ with the relative sign of the two terms changed. With this notation, the correlation 
function corresponding to the Xi interpolating field is 



G^(t,p;T) = i^V^-tr [-r [§ (s d , CS s C-\ S„) +g(s u , CS s C~\ S d ) } 



(24) 



and similarly for the correlation functions involving the X2 interpolating field. 



IV. LATTICE SIMULATIONS 

Having outlined the method of calculating excited 
baryon masses and the choice of interpolating fields, we 
next describe the gauge and fermion actions used in the 
present analysis. Additional details of the simulations 
can be found in Ref. |Tfij . 



A. Gauge Action 

For the gauge fields, the Luscher-Weisz mean-field im- 
proved plaquette plus rectangle action [3jJ is used. We 
define 

sq 

- i^£^ etr (l-EWz)), (25) 

rect 



where the operators U sq (x) and U YCC t(x) are defined as 

U sq (x) = U^,{x)U v {x + p)Ul{x + v)Ul{x) ,(26a) 
E/rcct(x) = U t _ l (x)U y (x + jx)U v (x + V + (j) 
xUl{x + 2v)Ut{x + v)Ul{x) 
+ U ft (x)U li (x + (i)U v (x + 2(i) 

xU£{x + ji + i>)Ut(x + 0)U}(x) .(26b) 

The link product U Iec t(x) denotes the rectangular 1x2 
and 2x1 plaquettes, and for the tadpole improvement 
factor we employ the plaquette measure 

/I \ V4 
u = ( -TZetT(U sq )) . (27) 

Gauge configurations are generated using the Cabibbo- 
Marinari pseudoheat-bath algorithm with three diagonal 
SU(2) subgroups looped over twice. Simulations are per- 
formed usin g a parallel algorithm with appropriate link 
partitioning |30|. 
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The calculations of octet excited-baryon masses are 
performed on a 16 3 x 32 lattice at {3 = 4.60. The scale is 
set via the string tension obtained from the static quark 
potential jU 



V(r) = V Q + . 



+ 1 



1 



where Vo, a, e and I are fit parameters, and [i] denotes 
the tree-level lattice Coulomb term 



r 






= 4 "/ 


r 





d 3 k 



•cos(k-r)L>44(0,k) 



with Z?44(fc) the time-time component of the gluon prop- 
agator. Note that Da{k^ k) is gauge-independent in the 
Breit frame, k± = 0, since k\jk 2 = 0. In the continuum 
limit. 



Taking the physical value of the string tension to be y/a 
440 MeV we find a lattice spacing of a = 0.122(2) fm. 



B. Fat-Link Irrelevant Fermion Action 

For the quark fields, we implement the Fat-Link Ir- 
relevant Clover (FLIC) action introduced in Ref. |16| . 
Fat links are created by averaging or smearing links on 
the lattice with their nearest transverse neighbours in a 
gauge covariant manner (APE smearing). In the FLIC 
action, this reduces the problem of exceptional configu- 
rations encountered with Wilson-style actions, and min- 
imises the effect of renormalisation on the action im- 
provement terms. By smearing only the irrelevant, higher 
dimensional terms in the action, and leaving the relevant 
dimension-four operators untouched, we retain short dis- 
tance quark an d gl uon interactions. Furthermore, the 
use of fat links 17] in the irrelevant operators removes 
the need to fine tune the clover coefficient in removing all 
0(a) artifacts. It is now clear that FLIC fermions pro- 
vide a new form of nonperturbative 0(a) improvement 
0- 

The smearing procedure |32( replaces a link, U^x), 
with a sum of the link and a times its staples 



U^x)^%{x) = {l-a)U^x) (28) 

4 

+ ^Y1 [ U A X ) U A X + va)Ul{x + fia) 

+Ul(x — ua)Un(x — ua)U v (x — va + fia) 

followed by projection back to SU(3). We select the uni- 
tary matrix [/J L which maximises 

KetriU^U'J) 



by iterating over the three diagonal SU(2) subgroups of 
SU(3). This procedure of smearing followed immediately 
by projection is repeated n times. The fat links used in 
this investigation are created with a = 0.7 and n = 4 
as discussed in Ref. 16] . The mean-field improved FLIC 
action is given by [1^] 



C FL e FL iCswKT - 

Z(u ) 



(29) 



where is constructed using fat links, and Uq is cal- 
culated via Eq. I|27|) using the fat links. The factor Csw 
is the (Sheikholeslami-Wohlert) clover coefficient [3j|, de- 
fined to be 1 at tree-level. The quark hopping parameter 
is k = l/(2m + 8r). We use the conventional choice of 
the Wilson parameter, r — 1. The mean-field improved 
Fat-Link Irrelevant Wilson action is 



Ufa - A) 



\ u 



uo 

u^Hx-n) 



rLr ' l,] ) ~ r ( ~ FL - ^(a + A) 



,FL 



(30) 



Our notation for the fermion action uses the Pauli rep- 
resentation of the Dirac 7-matrices defined in Appendix 
B of Sakurai |34j . In particular, the 7-matrices are Her- 
mitian with a^ v = h^, j v ]/(2i). 

As shown in Ref. |16j | , the mean-field improvement pa- 
rameter for the fat links is very close to 1, so that the 
mean-field improved coefficient for Csw is adequate ^(| . 
Another advantage is that one can now use highly im- 
proved definitions of F^ v (involving terms up to Uq 2 ), 
which give im pre ssive near-integer results for the topo- 
logical charge |35l | . 

In particular, we employ an 0(a ) improved definition 
of F„„ in which the standard clover-sum of four lxl 
loops lying in the fi, v plane is combined with 2x2 and 
3x3 loop clovers. Bilson- Thompson et al. [35| find 



p — 



-w 1 

2 

1 

90uf 



20u% 



(31) 



W 



3x3\ 



h.C 



tracclcss 



where W nxn is the clover-sum of four n x n loops and 
Fpu is made traceless by subtracting 1 /3 of the trace from 
each diagonal element of the 3x3 colour matrix. This 
definition reproduces the continuum limit with 0(a 6 ) er- 
rors. 

A fixed boundary condition in the time direction is 
used for the fermions by setting U t (x,N t ) = V x 
in the hopping terms of the fermion action, with peri- 
odic boundary conditions imposed in the spatial direc- 
tions. Gauge-invariant gaussian smearing 36] in the spa- 
tial dimensions is applied at the source to increase the 
overlap of the interpolating operators with the ground 



7 



states. The source-smearing technique [36j starts with a 
point source, ipo(xo,to), at space-time location (x ,£ ) = 
(1,1,1,3) and proceeds via the iterative scheme, 



i> i {x,t)=J2 F ( x > x ')^i-i( x '> t ), (32) 



where 



1 3 
F< ~ X ' X '' ) = (1 + a) V x ' x ' + f ^ [ U m( x ) 5 x',x+n 

* ' jX= 1 

+[/;(x-/2)^,^]) . (33) 



Repeating the procedure N times gives the resulting 
fermion field 

i> N (x,t) = Y,F N {x,x')M x ',t)- (34) 



The parameters N and a govern the size and shape of the 
smearing function and in our simulations we use N = 20 
and a = 6. 

Five masses are used in the calculations and the 
strange quark mass is taken to be the second heaviest 
quark mass in each case. The analysis is based on a 
sample of 400 configurations, and the error analysis is 
performed by a third-order, single-elimination jackknife, 
with the x 2 P er degree of freedom (TVdf) obtained via 
covariance matrix fits. 



V. CORRELATION MATRIX ANALYSIS 

In this section we outline the correlation matrix formal- 
ism for calculations of masses, coupling strengths and op- 
timal interpolating fields. After demonstrating that the 
correlation functions are real, we proceed to show how a 
matrix of such correlation functions may be used to iso- 
late states corresponding to different masses, and also to 
give information about the coupling of the operators to 
each of these states. 



A. The U + U* method 

A lattice QCD correlation function for the operator 
XiXjy where Xi is the i-th interpolating field for a par- 
ticular baryon (e.g. x% + m Section III), can be written 



g tj = (n|r(xiXi)|n) 



(35) 



XiXj 



where spinor indices and spatial coordinates are sup- 
pressed for ease of notation. The fermion and gauge ac- 
tions can be separated such that S[U, V», ip] = Sg[U] + 



i/jM[U]ip. Integration over the Grassmann variables ip 
and tp then gives 



/ VUe~ Sa M det(M[U])H t3 [U] 
JVUe-Selu] det(M[C/]) 



(36) 



where the term stands for the sum of all full contrac- 
tions of XiXj ■ The pure gauge action Sg and the fermion 
matrix M satisfy 



and 



S G [U] = S G [U*} 



TM[Z7*]r _1 = M*[U] 



(37) 



(38) 



respectively, where T is C75 in the Sakurai convention, 
adopted in Sec. IV addressing the lattice actions. 
Using the result of Eq. (|38|l . one has 



det(M [£/*]) = det(Af*[Z7]) 
and since det(M[E/]) is real, 

det (M[U*]) = det [M[U]) . 



(39) 



(40) 



Thus, U and U* are configurations of equal weight in the 
measure J X>f7det(M[£7])exp (— Sg[U]), in which case Qij 
can be written as 



On = 



1 / JVUe' s ^ det(M[J7]) {H i3 [U] + H^U*]} 



Let us define 



JVUe-Sciu] det(M[J7]) 



Gf 3 = tr sp {r ± g ij } , 



(41) 



(42) 



where tr sp denotes the spinor trace and r± is the 
parity-projection operator defined in Eq. Q. If 
tr sp {THijiU*}} = tr sp {TH* 3 [[/]}, then G% is real. This 
can be shown by first noting that ify will be products 
of 7-matrices, fermion propagators, and link-field oper- 
ators. Fermion propagators have the form M~ l and 
recalling that since TM[U*}T- 1 =M*[U], then we have 
TA/- 1 [[/*]r- 1 =(M" 1 [U])*. For products of link-field 
operators 0[U] contained in Hij, the condition 0[J7*] = 
0*[J7] is equivalent to the requirement that the coeffi- 
cients of all link-products are real. As long as this re- 
quirement is enforced, we can then simply proceed by 
inserting TT~ X inside the trace to show that the (spinor- 
traced) correlation functions are real. 

In summary, provided that the products of link op- 
erators in the interpolating fields are constructed using 
only real coefficients then the correlation functions Gfj 
are real. This symmetry is explicitly implemented by in- 
cluding both U and U* in the ensemble averaging used to 
construct the lattice correlation functions, providing an 
improved unbiased estimator which is strictly real. This 
is easily implemented at the correlation function level by 
observing 

m-\{u;}) = [C To M-\{U,}) (C75)- 1 ]* 

for quark propagators. 
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B. Recovering masses, couplings and optimal 
interpolators 

Let us again consider the momentum-space two-point 
function for t > 0, 

g l3 {t,p) = e-^ s (n\ Xi (t, <z)xj(o, o)|0) . (43) 

X 

At the hadronic level, 

x p' ,s B 

x (B,p', S |x,-(0,0)|fi) , 



It is convenient in the following discussion to label the 
states which have the x interpolating field quantum num- 
bers and which survive the parity projection as \B a ) for 
a = 1, 2, • • • , N. In general the number of states, N, in 
this tower of excited states may be infinite, but we will 
only ever need to consider a finite set of the lowest such 
states here. After selecting zero momentum, p = 0, the 
parity-projected trace of this object is then 

N 

G±(t) = tr sp {r ± g y (<, 0)} = ]T e-^Xfj; , (46) 

a=l 

— <y 

where Xf and A^ are coefficients denoting the couplings 
of the interpolating fields Xi an d Xj > respectively, to the 
state \B a ). If we use identical source and sink interpo- 
lating fields then it follows from the definition of the cou- 
pling strength that Xj = (Xf)* and from Eq. (|46JI we see 
that GfAt) = [G^ (*)]*, i.e., G ± is a Hermitian matrix. 
If, in addition, we use only real coefficients in the link 
products, then G ± is a real symmetric matrix. For the 
correlation matrices that we construct we have real link 
coefficients but we use smeared sources and point sinks 
and so in our calculations G is a real but non-symmetric 
matrix. Since G is a real matrix for the infinite number 
of possible choices of interpolating fields with real coeffi- 
cients, then we can take Xf and A^ to be real coefficents 
here without loss of generality. 

Suppose now that we have M creation and annihilation 
operators, where M < N. We can then form an M x M 
approximation to the full N x N matrix G. At this point 
there are two options for extracting masses. The first is 
the standard method for calculation of effective masses 



where the \B,p',s) are a complete set of states with mo- 
mentum p' and spin s 



EEE^y-^'^H 7 - (44) 

0> B s 



We can make use of translational invariance to write 



(45) 



at large t described in Section II. The second option is 
to extract the masses through a correlation-matrix pro- 
cedure |37| . 

Let us begin by considering the ideal case where we 
have N interpolating fields with the same quantum num- 
bers, but which give rise to N linearly independent states 
when acting on the vacuum. In this case we can constrict 
N ideal interpolating source and sink fields which per- 
fectly isolate the N individual baryon states \B a ), i.e., 



N 

8=1 

JV 

r = X, , (47b) 

such that 

{B P \T \Q) = 5 af3 z a u(a,p,s) , (48a) 
{Sl\<t> a \Bfi) = 5 aP z a u(a,p,s) , (48b) 

where z a and ~z a are the coupling strengths of 4> a and <fi 
to the state \B a ). The coefficients uf and v* a in Eqs. l|T7|) 
may differ when the source and sink have different smear- 
ing prescriptions, again indicated by the differentiation 
between z a and z a . For notational convenience for the 
remainder of this discussion repeated indices i,j,k are 
to be understood as being summed over. At p = 0, it 



e Ht e -iP.3 x .^ e iP.g e -m 



x f? s B » 

= E E e ~ EBt <nixi(o)|B,p, s)(B, P , S |xi(o)|o) 

i 



B,p',s) (B,p',s\ Xj (0)\n) 



s B 
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follows that, 



G % (*) u ? = (x>*p { r ± «? 



(49) 



The only i-dependence in this expression comes from the 
exponential term, which leads to the recurrence relation- 
ship 



which can be rewritten as 



[GHt + lTnGfAt) 



k J 



l k ■ 



(50) 



(51) 



This is recognized as an eigenvalue equation for the ma- 
trix [G ± (t + l)]- 1 G ± {t) with eigenvalues e m ° and eigen- 
vectors u a . Hence the natural logarithms of the eigenval- 
ues of [G ± (t+l)] _1 G ± (t) are the masses of the N baryons 
in the tower of excited states corresponding to the se- 
lected parity and the quantum numbers of the \ fields. 
The eigenvectors are the coefficients of the \ fields pro- 
viding the ideal linear combination for that state. Note 
that since here we use only real coefficients in our link 
products, then [G±(i + l)]- 1 G ± {t) is a real matrix and 
so u a and v a will be real eigenvectors. It also then fol- 
lows that z a and z" will be real. These coefficients are 
examined in detail in the following section. 

One can also construct the equivalent left-eigenvalue 
equation to recover the v vectors, providing the optimal 
linear combination of annihilation interpolators, 



vl a Gt j (t)=e m «vrGf j (t+l). 
Recalling Eq. (|49|l . one finds: 



(52) 



G?j(t) u« = r'Xfe-^ , (53) 



< a GfAt) = z a \j e~ 



(54) 



vT G± (t)G±(t) uf = z^A^V 2 "^ . (55) 



The definitions of Eqs. Q48|l imply 



(56) 



indicating the eigenvectors may be used to construct a 
correlation function in which a single state mass m a is 
isolated and which can be analysed using the methods of 
Section II. We refer to this as the projected correlation 
function in the following. Combining Eqs. I|55|) and l|56l) 
leads us to the result 



v* k a G k] (t)Gg(t) uf 
vl a G kl {t)uf 



= AfAje"™"* . (57) 



By extracting all N 2 such ratios, we can exactly recover 
all of the real couplings Af and A^ of x% an d Xj respec- 
tively to the state \B a ). Note that throughout this sec- 
tion no assumptions have been made about the symme- 
try properties of G^j- . This is essential due to our use of 
smeared sources and point sinks. 



In practice we will only have a relatively small num- 
ber, M < N, of interpolating fields in any given analysis. 
These M interpolators should be chosen to have good 
overlap with the lowest M excited states in the tower and 
we should attempt to study the ratios in Eq. (|57|l at early 
to intermediate Euclidean times, where the contribution 
of the (N — M) higher mass states will be suppressed but 
where there is still sufficient signal to allow the lowest M 
states to be seen. This procedure will lead to an esti- 
mate for the masses of each of the lowest M states in the 
tower of excited states. Of these M predicted masses, 
the highest will in general have the largest systematic 
error while the lower masses will be most reliably de- 
termined. Repeating the analysis with varying M and 
different combinations of interpolating fields will give an 
objective measure of the reliability of the extraction of 
these masses. 

In our case of a modest 2x2 correlation matrix 
(M = 2) we take a cautious approach to the selection 
of the eigenvalue analysis time. As already explained, we 
perform the eigenvalue analysis at an early to moderate 
Euclidean time where statistical noise is suppressed and 
yet contributions from at least the lowest two mass states 
is still present. One must exercise caution in performing 
the analysis at too early a time, as more than the desired 
M = 2 states may be contributing to the 2x2 matrix of 
correlation functions. 

We begin by projecting a particular parity, and then 
investigate the effective mass plots of the elements of the 
correlation matrix. Using the covariance-matrix based 
X 2 /-/Vdf, we identify the time slice at which all correla- 
tion functions of the correlation matrix are dominated by 
a single state. In practice, this time slice is determined 
by the correlator providing the lowest-lying effective mass 
plot. The eigenvalue analysis is performed at one time 
slice earlier, thus ensuring the presence of multiple states 
in the elements of the correlation function matrix, min- 
imising statistical uncertainties, and hopefully providing 
a clear signal for the analysis. In this approach mini- 
mal new information has been added, providing the best 
opportunity that the 2x2 correlation matrix is indeed 
dominated by 2 states. The left and right eigenvectors are 
determined and used to project correlation functions con- 
taining a single state from the correlation matrix as in- 
dicated in Eq. I|56|l . These correlation functions are then 
subjected to the same covariance-matrix based xV-^df 
analysis to identify new acceptable fit windows for deter- 
mining the masses of the resonances. 



VI. RESULTS 
A. Effective masses and the correlation matrix 

The correlation matrix analysis has a significant im- 
pact on the resolution of states obtained with the A c in- 
terpolating fields of Eqs. I|21|) and l|22|) . Hence we begin 
our discussion with a focus on these correlation functions. 
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FIG. 1: Effective masses of the lowest lying positive and neg- 
ative parity A states obtained using the A c interpolating field 
from 400 configurations using the FLIC action defined with 4 
sweeps of smearing at a = 0.7. The J p — i + (~ ) states la- 
beled Af (Af*) and A£ (A") are obtained using the XiXi an d 
X2X 2 interpolating fields, respectively. The smeared source is 
at t — 3. 



The effective mass plots for the positive and negative 
parity A states obtained using the A c interpolating field 
in the X1X1 and X2X2 correlation functions are shown in 
Fig. 1 for the FLIC action. Good values of the covariance 
matrix based xV-^df are obtained for the ground state 
(Af) for many different time-fitting intervals as long as 

one fits after time slice 9. Similarly, the lowest J p = | 
excitation for the X1X1 correlator (AJ*) requires fits fol- 
lowing time slice 8. The ground state (Af) mass ob- 
tained from X1X1 alone uses time slices 10-14 while the 
first odd-parity excited state, (Af*) uses time slices 9-12. 
The states obtained from the X2X2 correlation function 
plateau at earlier times and are also subject to noise ear- 
lier in time than the states obtained with the X1X1 cor- 
relator. For these reasons, good values of xV-^df are 
obtained on the time interval 6-8 for the positive par- 
ity states (A2), and time interval 8-11 for the negative 
parity states (Af). Hence, the time slice at which the 
eigenvalue analysis of the correlation matrix is performed 
is at T — 9 for the even-parity pair of states and at T = 8 
for the odd-parity pair of states. Selecting only one time 
slice earlier than that allowed by x 2 considerations pro- 
vides the best chance that only two states are present in 
the correlation matrix at that time. 

To guarantee the robustness of the eigenvector analysis 
and the subsequent projection procedure, various consis- 
tency checks are made at each stage of the process. For 



FIG. 2: As in Fig. 1, but for states obtained using the cor- 
relation functions projected from the correlation matrix as in 
Eq. ra . 



instance, a check is made to determine that the eigen- 
value in Eq. lj5U|) is positive, and that the mass deter- 
mined from the projected correlation function defined in 
Eq. 156( 1 is within the statistical fluctuations of the mass 
extracted without this analysis. For the octet interpolat- 
ing fields, off-diagonal elements are often suppressed by 
an order of magnitude relative to the diagonal elements 
and statistical noise can prevent the eigenvalue analysis 
from being successful. However, the strong suppression 
of off-diagonal elements is a clear signature that the mix- 
ing of the interpolating fields in these states is negligible. 

When the consistency checks are not satisfied, we have 
explored the possibility of stepping back to the previous 
time slice and performing the correlation matrix analysis 
there. In some cases, the mass of the lower-lying state 
reliably obtained via Euclidean time evolution is seen to 
increase in the eigenvalue analysis, indicating a failure 
of the correlation matrix analysis. The increase in the 
eigenvalue indicates that there are significant contribu- 
tions from three or more states in the 2x2 correlation 
matrix, thus spoiling the possibility of successful state 
isolation. In this case, the correlation matrix analysis is 
unable to provide additional information and masses are 
reported from the XiXi or X2X2 correlators as appropri- 
ate. 

Figure |21 illustrates the effective mass plots of the cor- 
relation functions projected from the correlation matrix 
as in Eq. (|56l) . The improved plateau behaviour is readily 
visible. Whereas in Fig. 1 the odd-parity effective masses 
are crossing at t = 6 and have minimal mass splitting, 
significant mass splitting between the two states is al- 
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TABLE I: Values of k used in this analysis and the corresponding pion and nucleon resonance masses for the FLIC action with 
4 sweeps of smearing at a = 0.7. Here K cr = 0.1300, and a string tension analysis gives a = 0.122(2) fm for ^/a — 440 MeV. 



mM 2 a 



0.1260 
0.1266 
0.1273 
0.1279 
0.1286 



0.5807(18) 
0.5343(19) 
0.4758(21) 
0.4203(23) 
0.3457(28) 



1.0972(49) 
1.0400(53) 
0.9701(59) 
0.9067(67) 
0.8273(86) 



1.388(14) 
1.340(16) 
1.286(19) 
1.244(25) 
1.186(33) 



1.442(12) 
1.404(15) 
1.363(20) 
1.345(29) 
1.374(57) 



1.676(12) 
1.642(13) 
1.605(15) 
1.580(18) 
1.571(26) 
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FIG. 3: As in Fig. 1, but for the nucleon states obtained using 
the correlation functions defined in Eqs. @ and 11UI . 



ready apparent at t — 6 in Fig. [2] The covariance based 
X 2 /A^df indicates that acceptable plateaus in the effec- 
tive mass plots start even earlier in some cases. The in- 
crease in mass splitting between the two negative parity 
states is more dramatic for A* than for the octet baryon 
interpolating fields. There the off diagonal elements of 
the correlation matrix are suppressed for the negative 
parity octet baryons, but not so for A*. As a result, 
the projection of states has only only a small effect for 
the octet baryon interpolators and this is detailed in Sec- 
tion Vic. 

Figs. |3] and 0] show the effective mass plots of the nu- 
cleon correlation functions XlXi an< l X2X2 an d following 
projection of the correlation matrix, respectively. Plots 
for the lightest quark mass considered are presented. The 
covariance matrix analysis of all quark masses indicates 
the following analysis windows in Euclidean time, 

N t , 10 - 14 ; N* , 9 — 12 ; 
A/J,8-ll ; N 2 ,6- 8. 



FIG. 4: As in Fig. 2, but for the nucleon states obtained 
using the correlation functions projected from the correlation 
matrix as in Eq. 156H . 



A comparison of Figs. El and 0| indicates that the cor- 
relation matrix analysis has a significantly smaller effect 
for the nucleon interpolators than the A c interpolators. 
This suggests that the states created by the interpolating 
fields Xi an d Xi have good overlap with the two lowest- 
lying physical nucleon states. 



B. Resonance masses and lattice action dependence 

In Fig. 5 we show the nucleon and N*(i ) masses 
as a function of the pseudoscalar meson mass squared, 
to 2 . The results of the new simulations are indicated by 
the filled squares for the FLIC action, and by the stars 
for the Wilson action (the Wilson points are obtained 
from a sample of 50 configurations). The values of m 2 
correspond to k values given in Table [I] 

For comparison, we also show results from earlier simu- 
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FIG. 5: Masses of the nucleon (N) and the lowest J p = k excitation ("JV*"). The FLIC and Wilson results are from the 
present analysis, with the DWF f2^l and NP improved clover |23| results shown for comparison. The empirical nucleon and 
low lying N*(^ ) masses are indicated by the asterisks along the ordinate. 



lations with domain wall fermions (DWF) [22| (open tri- 
angles), and a nonperturbatively (NP) improved clover 
action at /3 = 6.2 [23]. The scatter of the different NP 
improved results is due to different source smearing and 
volume effects: the open squares are obtained by using 
fuzzed sources and local sinks, the open circles use Ja- 
cobi smearing at both the source and sink, while the open 
diamonds, which extend to smaller quark masses, are ob- 
tained from a larger lattice (32 3 x 64) using Jacobi smear- 
ing. The empirical masses of the nucleon and the three 
lowest I excitations are indicated by the asterisks along 
the ordinate. 

There is excellent agreement between the different im- 
proved actions for the nucleon mass, in particular be- 
tween the FLIC, DWF ;22| and NP improved clover 
results. On the other hand, the Wilson results lie sys- 
tematically low in comparison to these due to the large 
0(a) errors in this action [16j. A similar pattern is re- 
peated for the N*{\~) masses. Namely, the FLIC, DWF 
and NP improved clover masses are in good agreement 
with each other, while the Wilson results again lie sys- 
tematically lower. A mass splitting of around 400 MeV 
is clearly visible between the N and N* for all actions, 
including the Wilson action, despite its poor chiral prop- 
erties. Furthermore, the trend of the iV*(| ) data with 



decreasing is consistent with the mass of the lowest 
lying physical negative parity N* states. 

Figure 6 shows the mass of the J p = \ + states (the ex- 
cited state is denoted by 'W(l/2 + )"). As is long known, 
the positive parity \i interpolating field does not have 
good overlap with the nucleon ground state and the 
correlation matrix results confirm this result, as discussed 
below. It has been speculated that \2 may have overlap 
with the lowest | + excited state, the TV* (1440) Roper 
resonance 22] . In addition to the FLIC and Wilson re- 
sults from the present analysis, we also show in Fig. 6 
the DWF results j2^| , and results from an earlier analysis 
with Wilson fermions together with the operator product 
expansion [l]|- The physical values of the lowest three 
i + excitations of the nucleon are indicated by the aster- 
isks. 

The most striking feature of the data is the relatively 
large excitation energy of the A'(i + ), some 1 GeV above 
the nucleon. There is little evidence, therefore, that this 
state is the TV* (1440) Roper resonance. While it is pos- 
sible that the Roper resonance may have a strong non- 
linear dependence on the quark mass at to 2 < 0.2 GeV 2 , 
arising from, for example, pion loop corrections, it is un- 
likely that this behaviour would be so dramatically dif- 
ferent from that of the TV* (1535) so as to reverse the level 
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ordering obtained from the lattice. A more likely expla- 
nation is that the \i interpolating field does not have 
good overlap with either the nucleon or the TV* (1440), 
but rather (a combination of) excited i + state(s). 

Recall that in a constituent quark model in a harmonic 
oscillator basis, the mass of the lowest mass state with 
the Roper quantum numbers is higher than the lowest 
P-wave excitation. It seems that neither the lattice data 
(at large quark masses and with our interpolating fields) 
nor the constituent quark model have good overlap with 
the Roper resonance. Better overlap with the Roper is 
likely to require more exotic interpolating fields. 

In Fig. 7 we show the ratio of the masses of the low- 
lying iV*(i ) and the nucleon. Once again, there is good 
agreement between the FLIC and DWF actions. How- 
ever, the results for the Wilson action lie above the oth- 
ers, as do those for the anisotropic D234 action [2l|. The 
D234 action has been mean-field improved, and uses an 
anisotropic lattice which is relatively coarse in the spatial 
direction (a « 0.24 fm). This is perhaps an indication of 
the need for nonperturbative or FLIC improvement. 



C. Resolving the resonances 

The mass splitting between the two lightest N*(h ) 
states (TV* (1535) and iV*(1650)) can be studied by con- 
sidering the odd parity content of the Xi an d Xi inter- 
polating fields in Eqs. JHJ and Q. Recall that the "di- 
quarks" in xx an d Xi couple differently to spin, so that 
even though the correlation functions built up from the 
Xi and X2 fields will be made up of a mixture of many 
excited states, th ey w ill have dominant overlap with dif- 
ferent states |2(|- By using the correlation-matrix 
techniques introduced in the previous section, we extract 
two separate mass states from the Xx an d X2 interpo- 
lating fields. The results from the correlation matrix 
analysis are shown by the filled symbols in Fig. 8 and 
are compared to the standard "naive" fits performed di- 
rectly on the diagonal correlation functions, xxXx an d 
X1X11 indicated by the open symbols. 

The results indicate that indeed the 7V*(i ) largely 
corresponding to the X2 field (labeled "JV|") lies above 
the N* ( I ) which can also be isolated via Euclidean time 
evolution with the xi field ("JV*") alone. The masses of 
the corresponding positive parity states, associated with 
the xi an d Xi fields (labeled U N±" and "-ZV2", respec- 
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tively) are shown for comparison. For reference, we also 
list the experimentally measured values of the low-lying 
I states. It is interesting to note that the mass splitting 
between the positive parity N\ and negative parity N* 2 
states (roughly 400-500 MeV) is similar to that between 
the N* 2 and the positive parity N2 state, reminiscent of 
a constituent quark-harmonic oscillator picture. 

The interpolating coefficients for the two positive and 
negative parity states (see Eq. 147(1 h extracted via the 
procedure outlined in Section IV Bl are given in Tables UTI 
and IIIII for various k values. The coefficients correspond- 
ing to each mass state (labeled "a" or "6" ) are normalised 
so that the sum of their absolute values is 1, 

KI + KI = 1 , (58a) 
KI + K4I = 1 , (58b) 

and similarly for the coefficients u"' 2 * for the negative 
parity mass states. This normalisation allows one to 
readily identify the fraction of each interpolating field 
needed to construct a linear combination having maxi- 
mum overlap with a particular baryon state. The last 
column in Tables UTI and IIIII shows the time slice T where 
the correlation matrix eigenvalue analysis is performed. 
From Table ITTI one immediately sees that the coefficient 



K2, reflecting the fraction of xi required to isolate the 
ground state nucleon, is extremely small. This further 
supports the earlier observation that the \i interpolat- 
ing field does not have good overlap with the nucleon 
ground state. Tabic IIIII shows the coefficients for iso- 
lating the two lowest-energy negative-parity N* states 
using the \i an d X2 interpolating fields. A significant 
amount of mixing is observed between the two interpo- 
lating fields for the lower energy state, particularly at 
heavy quark masses. This result is anticipated by the 
long Euclidean time evolution required to achieve an ac- 
ceptable xV-^df for the N* effective mass illustrated in 
Fig. The higher N 1 / 2 state, however, is dominated 
by the \2 field, thus explaining the good effective mass 
plateau observed in Fig.[3]without the correlation matrix 
approach. Note that the most significant contribution to 
the N% state from xi is f° r the third quark mass when 
the correlation matrix analysis is performed at an early 
time slice and is spoiled by contamination from higher 
excited states. The most significant contribution at the 
preferred time slice, which also has the smallest errors, is 
for the lightest quark mass. It is for these reasons that 
we choose the lightest quark mass in Fig. [3] to illustrate 
the effective masses of the projected nucleon states. 
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FIG. 8: Masses of the J p = | and | nucleon states, for the FLIC action. The positive (negative) parity states are labeled 
Ni (N*) and N2 (iV|). The results from the projection of the correlation matrix as discussed in Sec. IVI K\ are shown by the 
filled symbols, whereas the results from the standard fits to the X1X1 an d X2X2 correlation functions are shown by the open 
symbols (offset to the right for clarity). Empirical masses of the low lying | states are indicated by the asterisks. 



Turning to the strange sector, in Fig. 9 we show the 
masses of the positive and negative parity £ baryons cal- 
culated from the FLIC action compared with the physical 
masses of the known positive and negative parity states. 
The data for the masses of these states are listed in Ta- 
ble ^3 and the interpolator coefficients for the two posi- 
tive and negative parity states are given in Tables |V| and 
IVII respectively. The pattern of mass splittings is simi- 
lar to that found in Fig. 8 for the nucleon. Namely, the 
i + state associated with the xi field appears consistent 

with the empirical £(1193) ground state, while the \ + 
state associated with the \i field lies significantly above 
the observed first (Roper-like) i + excitation, £*(1660). 
There is also evidence for a mass splitting between the 
two negative parity states, similar to that in the non- 
strange sector. The behaviour of the interpolator coef- 
ficients for the I] 1 / 2 and E 1 / 2 states is also similar to 
that for the nucleon in Tables |H] and 11111 Namely, while 
the positive parity ground state is dominated by the xi 
interpolating field, there is considerable mixing between 
the Xi an d X2 fields for the lowest negative parity state, 
with the higher E 1 / 2 state receiving a dominant contri- 
bution from xi- 

The spectrum of the strangeness -2 positive and nega- 



tive parity 3 hyperons is displayed in Fig. 10, with data 
given in Table IVIII and the interpolator coefficients for 
the ~}/ 2+ and S 1 / 2 " states in Tables IVTnl and El re- 
spectively. Once again, the pattern of calculated masses 
repeats that found for the £ and N masses in Figs. 8 and 
9, and for the respective coupling coefficients. The em- 
pirical masses of the physical 2* baryons are denoted by 
asterisks. However, for all but the ground state S(1318), 
the J p values are not known. 

Finally, we consider the A hyperons. In Figs. UTI and 
1121 we compare results obtained from the A 8 and A c in- 
terpolating fields, respectively, using the two different 
techniques for extracting masses. The data are given in 
Tables DO and IXllTl respectively. A direct comparison be- 
tween the positive and negative parity masses for the A 8 
(open symbols) and A c (filled symbols) states extracted 
from the correlation matrix analysis, is shown in Fig. 13. 
A similar pattern of mass splittings to that for the N* 
spectrum of Fig. 8 is observed. In particular, the neg- 
ative parity A* state (diamonds) lies ~ 400 MeV above 
the positive parity Ai ground state (circles), for both the 
A 8 and A c fields. There is also clear evidence of a mass 
splitting between the A* (diamonds) and Aj (squares). 

Using the naive fitting scheme (open symbols in 
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FIG. 9: As in Fig. 8 but for the E baryons. 



TABLE II: Interpolating field coefficients for the two positive 
parity A 1 / 2 states. The time slice T at which the correlation 
matrix analysis is performed is indicated in the last column. 



TABLE IV: E baryon resonance masses. 



K 


uf 


U a 2 


u! 


u§ 


T 


0.1260 


0.999(1) 


0.001(1) 


0.154(28) 


0.846(28) 


7 


0.1266 


0.997(2) 


0.003(2) 


0.112(59) 


0.888(59) 


8 


0.1273 


0.996(2) 


0.004(2) 


0.083(74) 


0.917(74) 


8 


0.1279 


0.993(3) 


0.007(3) 


0.049(99) 


0.951(99) 


8 


0.1286 


0.989(3) 


0.011(3) 


0.066(81) 


0.934(81) 


7 



k mg 1 a mgj a 

0.1260 1.0765(50) 1.371(15) 
0.1266 1.0400(53) 1.340(16) 
0.1273 0.9966(57) 1.307(17) 
0.1279 0.9589(62) 1.281(20) 
0.1286 0.9149(72) 1.265(29) 



1.432(13) 1.665(12) 

1.404(15) 1.642(13) 

1.371(17) 1.617(14) 

1.349(21) 1.597(16) 

1.332(28) 1.580(19) 



TABLE III: Interpolating field coefficients for the two nega- 
tive parity A 1//2 states. The eigenvalues of the correlation 
matrix analysis indicate that excited states spoil the eigen- 
state isolation for k values 0.1273 and 0.1279. 



K 


uT 


uT 


u\* 


6* 
u 2 


T 


0.1260 


0.45(7) 


0.55(7) 


0.16(15) 


-0.84(15) 


8 


0.1266 


0.50(7) 


0.50(7) 


0.08(14) 


-0.92(14) 


8 


0.1273 


0.42(6) 


0.58(6) 


0.26(15) 


-0.74(15) 


7 


0.1279 


0.53(10) 


0.47(10) 


0.09(15) 


-0.91(15) 


7 


0.1286 


0.77(15) 


0.23(15) 


0.20(5) 


0.80(5) 


8 



Figs. 111! and 112 . misses the mass splitting between A* 



and Aj for the "common" interpolating field. Only af- 
ter performing the correlation matrix analysis is it possi- 
ble to resolve two separate mass states, as seen by the 
filled symbols in Fig. This may be an indication 
that the physics responsible for the mass splitting be- 
tween the negative parity A*(1670) and A*(1800) states 
is suppressed in the A c interpolating field. This is also 
evidenced by comparing the interpolator coefficients for 
the positive and negative parity A 8 and A c states, in Ta- 
bles EH and EH and lxTvl andlXVl respectively. While 
the couplings for the A 8 for both the positive parity states 
are similar to those for the nucleon and other hyperons, 
there is more prominent mixing for the case of the A c . In 
particular, there is notably stronger mixing for the higher 
mass negative parity state in the case of the A c compared 
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FIG. 10: As in Fig. 8 but for the S baryons. The J p values of the excited states marked with "?" are undetermined. 



TABLE V: Interpolating field coefficients for the two E x / 2+ TABLE VII: 5 baryon resonance masses. 



states. 



K 


uf 


u% 


u\ 


u\ 


T 


K 

0.1260 


ms, a 
1.0612(52) 


1.358(15) 


1.414(13) 


ms 2 a 
1.653(12) 


0.1260 


0.997(1) 


0.003(1) 


0.127(53) 


0.873(53) 


8 


0.1266 


1.0400(53) 


1.340(16) 


1.404(15) 


1.642(13) 


0.1266 


0.997(2) 


0.003(2) 


0.112(59) 


0.888(59) 


8 


0.1273 


1.0145(54) 


1.320(16) 


1.392(17) 


1.630(14) 


0.1273 


0.998(2) 


0.002(2) 


0.133(35) 


0.867(35) 


7 


0.1279 


0.9919(56) 


1.302(18) 


1.389(20) 


1.622(15) 


0.1279 


0.997(2) 


0.003(2) 


0.121(41) 


0.879(41) 


7 


0.1286 


0.9649(60) 


1.281(20) 


1.399(27) 


1.618(24) 


0.1286 


0.996(3) 


0.004(3) 


0.100(52) 


0.900(52) 


7 













TABLE VI: Interpolating field coefficients for the two E 1/2 
states. The eigenvalues of the correlation matrix analysis in- 
dicate that excited states spoil the eigenstate isolation for k 
values 0.1273 through 0.1286. 



K 


uf* 


uT 


u\* 


uf 


T 


0.1260 


0.47(7) 


0.53(7) 


0.11(13) 


-0.89(13) 


8 


0.1266 


0.50(7) 


0.50(7) 


0.08(14) 


-0.92(14) 


8 


0.1273 


0.38(5) 


0.62(5) 


0.35(14) 


-0.65(14) 


7 


0.1279 


0.42(7) 


0.58(7) 


0.30(17) 


-0.70(17) 


7 


0.1286 


0.52(13) 


0.48(13) 


0.17(22) 


-0.83(22) 


7 



with the corresponding A 8 state. The \2 contributes 



~ 80-90% of the strength compared to - 50-60% for 
the X2 ■ The interpolator coefficients are precisely deter- 
mined in the A c correlation matrix analysis. As for the 
other baryons, there is little evidence that the A2 (tri- 
angles) has any significant overlap with the first positive 
parity excited state, A* (1600) (c.f. the Roper resonance, 
Ar* (1440), in Fig.©. 

While it seems plausible that nonanalyticities in a chi- 
ral extrapolation |jj of N\ and N* results could eventu- 
ally lead to agreement with experiment, the situation for 
the A* (1405) is not as compelling. Whereas a 150 MeV 
pion- induced self energy is required for the N± , JV-* and 
Ai , 400 MeV is required to approach the empirical mass 
of the A*(1405). This may not be surprising for the octet 
fields, as the A* (1405), being an SU(3) flavour singlet, 
may not couple strongly to an SU(3) octet interpolating 
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FIG. 11: As in Fig. 8 but for the A states obtained using the A 8 interpolating field. 



TABLE VIII: Interpolating field coefficients for the two 3 
states. 



TABLE X: A baryon resonance masses from the octet, A 
interpolating field. 



K 


uf 


u% 


u\ 


u\ 


T 


K 




m A *a 


m A *a 




0.1260 


0.999(1) 


0.001(1) 


0.146(30) 


0.854(30) 


7 


0.1260 


1.0801(50) 


1.374(15) 


1.427(13) 


1.665(12) 


0.1266 


0.997(2) 


0.003(2) 


0.112(59) 


0.888(59) 


8 


0.1266 


1.0400(53) 


1.340(16) 


1.404(15) 


1.642(13) 


0.1273 


0.996(2) 


0.004(2) 


0.105(64) 


0.895(64) 


8 


0.1273 


0.9910(56) 


1.302(17) 


1.380(19) 


1.618(15) 


0.1279 


0.995(2) 


0.005(2) 


0.097(70) 


0.903(70) 


8 


0.1279 


0.9464(61) 


1.269(21) 


1.373(26) 


1.603(17) 


0.1286 


0.993(2) 


0.007(2) 


0.076(83) 


0.924(83) 


8 


0.1286 


0.8904(72) 


1.233(28) 


1.410(47) 


1.599(21) 



TABLE IX: Interpolating field coefficients for the two S 1 ^ 2 
states. The eigenvalues of the correlation matrix analysis in- 
dicate that excited states spoil the eigenstate isolation for k 
values 0.1273 through 0.1286. 



« 


uT 


uT 


u\* 


uf 


T 


0.1260 


0.48(8) 


0.52(8) 


0.13(16) 


-0.87(16) 


8 


0.1266 


0.50(7) 


0.50(7) 


0.08(14) 


-0.92(14) 


8 


0.1273 


0.38(5) 


0.62(5) 


0.32(13) 


-0.68(13) 


7 


0.1279 


0.42(6) 


0.58(6) 


0.22(13) 


-0.78(13) 


7 


0.1286 


0.49(7) 


0.51(7) 


0.09(11) 


-0.91(11) 


7 



field. Indeed, there is some evidence of this in Fig. 13. 



This large discrepancy of 400 MeV suggests that relevant 
physics giving rise to a light A* (1405) may be absent 
from simulations in the quenched approximation. The 
behaviour of the A^ 2 states may be modified at small 
values of the quark mass through nonlinear effects asso- 
ciated with Goldstone boson loops including the strong 
coupling of the A* (1405) to S7r and KN channels. While 
some of this coupling will survive in the quenched approx- 
imation, generally the couplings are modified and sup- 
pressed 0, |3g ■ It is also interesting to note that the K\ 
and A 2 masses display a similar behaviour to that seen 
for the S* and 'S, states, which are dominated by the 
heavier strange quark. Alternatively, the study of more 
exotic interpolating fields may indicate the the A* (1405) 
does not couple strongly to \i or Xi- Investigations at 
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FIG. 12: As in Fig. 8 but for the A states obtained using the A c interpolating field. 



TABLE XI: Interpolating field coefficients for the two positive 
parity A 8 states. 



K 


ut 


u% 


ui 


u\ 


T 


0.1260 


0.999(1) 


0.001(1) 


0.149(29) 


0.851(29) 


7 


0.1266 


0.997(2) 


0.003(2) 


0.112(59) 


0.888(59) 


8 


0.1273 


0.995(2) 


0.005(2) 


0.095(69) 


0.905(69) 


8 


0.1279 


0.993(2) 


0.007(2) 


0.070(85) 


0.930(85) 


8 


0.1286 


0.990(2) 


0.010(2) 


0.081(63) 


0.919(63) 


7 



TABLE XIII: A baryon resonance masses from the "com- 
mon", A c , interpolating field. 



K 


mA 1 a 


m A *a 


m A *a 


mA 2 a 


0.1260 


1.0815(50) 


1.334(13) 


1.408(12) 


1.662(11) 


0.1266 


1.0413(52) 


1.301(14) 


1.382(13) 


1.638(12) 


0.1273 


0.9920(56) 


1.262(16) 


1.356(16) 


1.611(12) 


0.1279 


0.9473(61) 


1.226(18) 


1.342(21) 


1.590(13) 


0.1286 


0.8912(73) 


1.181(21) 


1.357(33) 


1.570(15) 



TABLE XII: Interpolating field coefficients for the two nega- 
tive parity A 8 states. The eigenvalues of the correlation ma- 
trix analysis indicate that excited states spoil the eigenstate 
isolation for k values 0.1273 through 0.1286. 





a* 


a* 


u\* 


6* 


T 


0.1260 


0.46(8) 


0.54(8) 


0.16(16) 


-0.84(16) 


8 


0.1266 


0.50(7) 


0.50(7) 


0.08(14) 


-0.92(14) 


8 


0.1273 


0.40(6) 


0.60(6) 


0.27(14) 


-0.73(13) 


7 


0.1279 


0.49(8) 


0.51(8) 


0.12(13) 


-0.88(13) 


7 


0.1286 


0.47(8) 


0.53(8) 


0.19(13) 


-0.81(13) 


6 



lighter quark masses involving quenched chiral perturba- 
tion theory will assist in resolving these issues. 



VII. CONCLUSION 

We have presented the first results for the excited 
baryon spectrum from lattice QCD using an 0(a 2 ) im- 
proved Luscher-Weise gauge action and an 0(a)- 
improved Fat-Link Irrelevant Clover (FLIC) quark ac- 
tion in which only the links of the irrelevant dimension 
five operators are smeared |lfi[ . The FLIC action pro- 
vides a new form of nonperturbative 0(a) improvement 
in which 0(a) errors are eliminated and 0(a 2 ) errors are 
very small [4fjj| . The simulations have been performed on 
a 16 3 x 32 lattice at j3 = 4.60, providing a lattice spacing 
of a = 0.122(2) fm. The analysis is based on a set of 400 
configurations generated on the Orion supercomputer at 
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FIG. 13: Masses of the positive and negative parity A states, for the octet A 8 (open symbols) and "common" A c (filled symbols) 
interpolating fields with the FLIC action. The positive (negative) parity states labeled Ai (A*) and A2 (A|) are the two states 
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TABLE XIV: Interpolating field coefficients for the two posi- 
tive parity A c states. 





u? 


u a 2 


u\ 


u\ 


T 


0.1260 


1.000(2) 


0.000(2) 


0.282(51) 


-0.718(51) 


9 


0.1266 


0.997(2) 


0.003(2) 


0.291(55) 


-0.709(55) 


9 


0.1273 


0.994(2) 


0.006(2) 


0.278(26) 


-0.722(26) 


8 


0.1279 


0.990(2) 


0.010(2) 


0.279(18) 


-0.721(18) 


7 


0.1286 


0.983(3) 


0.017(3) 


0.278(13) 


-0.722(13) 


6 


TABLE XV: Interpolating field coefficients for the two nega- 


tive parity A c states 


The correlation matrix analysis is 


suc- 


cessful for all k values. 










uT 


«r 


u\* 


u\* 


T 


0.1260 


0.54(2) 


-0.46(2) 


0.23(3) 


0.77(3) 


8 


0.1266 


0.53(2) 


-0.47(2) 


0.27(3) 


0.73(3) 


8 


0.1273 


0.52(1) 


-0.48(1) 


0.33(3) 


0.67(3) 


8 


0.1279 


0.51(1) 


-0.49(1) 


0.39(3) 


0.61(3) 


8 


0.1286 


0.49(1) 


-0.51(1) 


0.47(4) 


0.53(4) 


8 



the University of Adelaide. 

Good agreement is obtained between the FLIC and 
other improved actions, including the nonperturbatively 
improved clover |2^| and domain wall fermion (DWF) 
[22J actions, for the nucleon and its chiral partner, with a 
mass splitting of ~ 400 MeV. Our results for the iV*(|~) 
improve on those using the D234 |2l| and Wilson actions. 
Despite strong chiral symmetry breaking, the results with 
the Wilson action are still able to resolve the splitting 
between the chiral partners of the nucleon. Using the 
two standard nucleon interpolating fields, we also confirm 
earlier observations 20] of a mass splitting between the 
two nearby i states. We find no evidence of overlap 
with the | + Roper resonance. 

In the strange sector, we have investigated the overlap 
of various A interpolating fields with the low-lying i 
states. Once again a clear mass splitting of ~ 400 MeV 
between the octet A and its parity partner is seen, with 
evidence of a mass splitting between the two low-lying 
odd-parity states We find no evidence of strong overlap 
with the ^ + "Roper" excitation, A*(1600). The empirical 
mass suppression of the A* (1405) is not evident in these 
quenched QCD simulations, possibly suggesting an im- 
portant role for the meson cloud of the A*(1405) and/or 
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a need for more exotic interpolating fields. 

We have not attempted to extrapolate the lattice re- 
sults to the physical region of light quarks, since the non- 
analytic behaviour of iV*'s near the chiral limit is not as 
well studied as that of the nucleon pf |8|, |4l| . It is vital 
that future lattice N* simulations push closer towards 
the chiral limit. On a promising note, our simulations 
with the 4 sweep FLIC action are able to reach relatively 
low quark masses (m q ~ 60-70 MeV) already. Our dis- 
cussion of quenching effects is limited to a qualitative 
level until the formulation of quenched chiral perturba- 
tion theory for | baryon resonances is established 01 
or dynamical fermion simulations are completed. Expe- 
rience suggests that dynamical fermion results will be 
shifted down in mass relative to quenched results, with 
increased downward curvature near the chiral limit @. 
It will be fascinating to confront this physics with both 
numerical simulation and chiral nonanalytic approaches. 

In order to further explore the origin of the Roper res- 
onances or the A* (1405), more exotic interpolating fields 
involving higher Fock states, or nonlocal operators should 



be investigated. Finally, the present N* mass analysis 
will be extended in future to include TV — > N* transi- 
tion form factors through the calculation of three-point 
correlation functions. 
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